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ABSTRACT 



Aims. We examine an Infrared Dark Cloud at high spatial resolution as a means to study rotation, outflow, and infall at the onset of 
massive star formation. 

Methods. The Infrared Dark Cloud (IRDC) 18223-3 was observed at 1.1 mm and 1.3 mm with the Submillimeter Array (SMA) and 
follow-up short spacing information was obtained with the IRAM 30 m telescope. Additional data were taken at 3 mm with the IRAM 
Plateau de Bure Interferometer (PdBI). 

Results. Submillimeter Array observations combined with IRAM 30 meter data in 12 CO(2-l) reveal the outflow orientation in the 
IRDC 18223-3 region, and PdBI 3 mm observations confirm this orientation in other molecular species. The implication of the 
outflow's presence is that an accretion disk is feeding it, so using line data for high density tracers such as C 18 0, N9H + , and CH3OH, 
we looked for indications of a velocity gradient perpendicular to the outflow direction. Surprisingly, this gradient turns out to be most 
apparent in CH3OH. The large size (28,000 AU) of the flattened rotating object detected indicates that this velocity gradient cannot be 
due solely to a disk, but rather from inward spiraling gas within which a Keplerian disk likely exists. The rotational signatures can be 
modeled via rotationally infalling gas. From the outflow parameters, we derive properties of the source such as an outflow dynamical 
age of ~37,000 years, outflow mass of ~13 M Q , and outflow energy of ~1.7 x 10 46 erg. While the outflow mass and energy are clearly 
consistent with a high-mass star forming region, the outflow dynamical age indicates a slightly more evolved evolutionary stage than 
previous spectral energy distribution (SED) modeling indicates. 

Conclusions. The orientation of the molecular outflow associated with IRDC 18223-3 is in the northwest-southeast direction and 
velocity gradients orthogonal to the outflow reveal a large rotating structure likely harboring an accretion disk within. We also present 
a model of the observed methanol velocity gradient. The calculated outflow properties reveal that this is truly a massive star in the 
making. These data present evidence for one of the youngest known outflow/infall/disk systems in massive star formation. A tentative 
evolutionary picture for massive disks is discussed. 

Key words, stars: formation - stars: individual (IRDC 18223-3) - stars: early type 



1. Introduction 

Massive stars (M>8M ) are paramount constituents of galax- 
ies, yet their formation is significantly less well understood 
than low-mass star formation. Theories of massive star for- 
mation predict massive accretion disks and indirect evi- 
dence for disks is prevalent, but observational evidence for 
such structures is rare. Outflows appear to be widespread in 



massive star formation regions (Shepherd & Churchwell 1996 
Zhang et al. 200 1| 2005; IBeufher et al. 2002bl) . and must be 
powered and fed by a (massive) accretion disk. Evidence 
of rotation and infall that hint at the presence of disks has 
been detected in several systems, although further study is re- 
quired (see reviews by IBeufher et al. 2007UCeseroni et al. 20071 
IZinnecker & Yorke 20071 . 

Several factors complicate observational disk studies of mas- 
sive star formation regions. Compared to low mass star for- 
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* The FITS files of the integrated continuum and line maps, as well 
as the first and second moment maps are available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat 7J/A+ A/ 



mation regions, they are much rarer and further away on av- 
erage. Additionally, massive stars evolve quickly and typically 
form deeply embedded within a core at the center of a cluster, 
thus obscuring the process from optical wavelengths and mak- 
ing it difficult to distinguish one protostar from another in such a 
busy region of star formation typically at kpc distances. Due to 
these complications, it is extremely difficult to directly observe 
massive accretion disks at sufficiently high spatial resolution to 
definitively identify the disk component. Nevertheless, today's 
interferometers are conducive to studies that probe size-scales 
comparable to the sizes we might expect for accretion disks 
around massive stars, and studies of disks in massive star for- 
mation regions have recently increased significantly in number 
(ICeseroni et al. 20071 ). 

The detection of circumstellar disks would provide strong 
support for the theory that massive star formation is simi- 
lar to a scaled-up version of low-mass star formation (e.g. 
lYorke & Sonnhalter 20021 IKrumholz et al. 20071 2009). These 
rotating disks, if present, would exhibit a velocity gradient per- 
pendicular to the established outflow direction. 

While previous disk studies in high-mass star forma- 
tion mainly concentrated on relatively evolved regions such 
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as high-mass protostellar objects (e.g. IRAS 20126+4104, 
ICesaroni et al. 19971 2005; G 192.1 6-3 .82,|Shepherd et al. 2002 
IRAS 18089-1732. IBeuther& Walsh 20081) . here we focus on 
one of the earliest evolutionary stages, namely an Infrared Dark 
Cloud. 

A survey of high-mass protostellar objects (HMPOs) as- 
sociated with Infrared Astronomical Satellite (IRAS) sources 
was conducted (IBeuther et al. 2002al l. and IRDC 18223-3 hap- 
pened to be in the field of view of one of the target sources, 
IRAS 18223-1243. IR DC 18223-3 lies at a distance of 3.7 kpc 
( Sridharan et al. 20051 ) in a filamentary structure extending south 
from IRAS 18223-1243 ( IBeuther et al. 2005bl l. IRDC 18223-3 
is likely a high-mass core containing a young stellar object that 
is currently accreting at a high enough rate that it will even- 
tually join the high-mass regime (Beuther & Steinacker 2007). 
Our current studies present the outflow properties of this source 
and reveal velocity gradients in several molecular species that 
are not aligned with the outflow orientation. 

Based on the analysis of Spitzer observations combined 
with SMA observations, IRDC 18223-3 is comprised of at 
least two separate components (Beuther & Steinacker 2007]). A 
cooler component at ~15 K contains most of the mass (580 M ) 
while a warmer component at ~50 K with approximately 0.01 
M is necessary to explain the excess flux in the SED at 24 fim. 
The young source has an accretion-dominated luminosity of 180 
L Q that is expected to increase over time. While the young stellar 
object has not yet accumulated a significant amount of mass, in 
the framework of the simulations by Krumholz et al. (2007), the 
source's mass and luminosity are consistent with a large accre- 
tion rate of ~10~ 4 M Q yr _1 . This high accretion rate, as well as 
the outflow parameters we derive, both point to a scenario that 
this source will eventually become a massive star and that we 
are looking at one of the earliest stages of massive star forma- 
tion. Therefore, this source is an excellent case study for looking 
at the properties of massive stars at the beginning stages of their 
lifetimes. This source also provides an ideal specimen to study 
the role of disks in early evolutionary phases of massive star for- 
mation. 

In this paper we report a cone-shaped structure detected in 
12 CO(2-l) outlining the molecular outflow as well as a velocity 
gradient of 3 km s' 1 in CH 3 OH(6^-5^) and N 2 H + (3-2) perpen- 
dicular to the outflow. We also see velocity gradients on the order 
of a few km s _I not aligned with the outflow axis in C 18 0(2-l). 
Additionally, we present the results of modeling the observed 
CH3OH parameters to obtain the velocity and density structure 
of the region. 

2. Observations 

Observations with the Submillimeter Array, Plateau de Bure 
Interferometer and IRAM 30 m observatories include the 1 .3 mm 
dust continuum and several spectral lines (see TableQ]). The de- 
tails of the observations are discussed below. 

2. 1 . Submillimeter Array 

Observations of IRDC 18223-3 with the Submillimeter AiTa>Q 
were carried out the nights of 30 May 2006, 8 August 2006, 24 
July 2007, and 27 August 2007. The phase reference center of 

1 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 



Table 1. Observed lines 
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this source is RA. 18h25m08.55s Decl.= -12°45'23.3" (J2000.0 ) 
and the velocity of rest v/ jr is 45.3 km s _1 (ISridharan et al. 20051 ). 
The SMA has two spectral sidebands-both 2 GHz wide sepa- 
rated by 10 GHz. All SMA data sets were calibrated with the 
IDL superset MIRQ, then imaged and analyzed with MIRIAD. 
These submillimeter interferometric observations allow us to 
probe through the optically opaque molecular cloud which IRDC 
18223-3 is deeply embedded within. Specific details of the ob- 
servations are listed in Table [2] 

The observations on 30 May 2006 are with all 8 antennas at 
1.3 mm in the extended configuration. Zenith opacity measured 
by the Caltech Submillimeter Observatory was approximately 
t(230) 0.15. Observing with 8 antennas coupled with the fact 
that the track was 8.5 hours long (5 hours on-source integration) 
provided optimal coverage in the uv plane. Passband and flux 
calibrations were derived from observations of Uranus. For the 
flux, solutions could not be found on all baselines so it was nec- 
essary to perform an initial flux calibration on 3C279 to set the 
relative flux levels. Phase and amplitude calibrations were made 
from regular observations (5 minutes every half hour) of quasars 
1911-201 and 1733-130. The same spectral setup was observed 
on 8 August 2006 in the compact configuration. Only 7 anten- 
nas were used in this data set. Despite t varying between 0.2 
and 0.25, the data are good quality. The same calibration objects 
were used as in the 30 May data set, however, it was not nec- 
essary to set the relative flux levels before using Uranus for the 
flux calibration. The track was 7 hours long, 4 hours of which 
were spent on-source. After combination of the data from both 
configurations, projected baselines range from 9 to 163 kA. 

Observations at 1.1 mm were obtained on 24 July 2007 in 
the compact configuration with all 8 antennas and in the ex- 
tended configuration with only 6 fully-functioning antennas on 
27 August 2007. Zenith opacity, t(230 GHz), was -0.15 for the 
compact configuration observations, and ~0. 17 for the extended 
configuration observations. The quasar 1911-201 was used for 
phase and amplitude calibrations for the compact configuration 
data, and quasar 1733-130 was used for the extended configura- 
tion data. Uranus was used for flux calibration in both 1 . 1 mm 
data sets and for phase and amplitude calibration of the compact 
configuration. Quasar 3C454.3 was used for bandpass calibra- 
tion of the extended configuration data as well as to set relative 
flux levels prior to using Uranus for the absolute flux calibration. 
Projected baselines ranged from 11 to 189 kA in the combined 
1 . 1 mm data set. 

2 The MIR cookbook by Charlie Qi can be found at 
http://cfa -www.harvard.edu/~cqi/mircook.html 
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Entries include the observatory at which the data were taken, the date of observation, the frequency observed, the array configuration, the 
calibrators, the synthesized beam obtained after inverting and cleaning the data, the l<x continuum noise level, and the 1 km s -1 rms. 
"Units: K at 0.4 km s -1 resolution averaged over the spectra in the central 4"x4" of the map. 



2.2. Plateau de Bure Interferometer 

Additional observations of IRDC 1 8223-3 were obtained at 3mm 
with the IRAM Plateau de Bure Interferometer (PdBlj^ located 
near Grenoble, France. Observations were carried out on 23 May 
2006, then 3 and 5 April 2007. These observations include de- 
tections of SiO(2-l), S0 2 , HN 13 C, and H 13 CO + (1-0) as well as 
the 3 mm continuum. Quasars NRAO 530 and 1741-038 were 
used as phase and amplitude calibrators and MWC 349 for flux 
calibration. Further details are listed in Table [2] 

2.3. Pico Veleta 30 m Telescope 

To complement our interferometric data, we obtained short spac- 
ing observations of IRDC 18223-3 on 24 February 2007 with 
the HERA receiver at the IRAM 30 m telescope near Granada, 
Spain. The 12 CO(2-l) line at 230.5 GHz and 13 CO(2-l) at 220.4 
GHz were observed in on-the-fly mode, resulting in maps ap- 
proximately 140" x 140" in extent. The sampling interval was 
2" and the region was scanned four times in the north-south di- 
rection and four times in the east-west direction in order to re- 
duce effects caused by the scanning process. The spectra were 
reduced using CLASS which is part of the GILDAS software 
package and then right ascension and declination scans were 
combined with the plait algorithm in GREG, another component 
of the GILDAS collection. The 12 CO data have a beam size of 
11" and a T m /, corrected rms noise level of 0. 17 K at 0.4 km s _1 
resolution (averaged over the spectra in the central 4"x4" region 
of the map). 

After reducing the 30 m data set separately, single dish 12 CO 
data were converted to visibilities and then combined with the 
SMA data sets using the MIRIAD package task UVMODEL. 
The synthesized beam of the combined data is 1.9" x 1.7" (PA 
70°). 



3. Results and Discussion 

3. 1 . Large Scale Emission 

IRDC 18223-3 is about 3.5' south of IRAS 18223-1243 in a fila- 
mentary structure extending more than 5' (>5 pc) roughly in the 
north-south direction across the sky. In the large scale map pre- 
sented in Figure [TJ the 1.2 mm continuum from previous obser- 
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Fig. 1. Image of the velocity integrated 12 CO(2-l) flux from the 
IRAM 30m single dish telescope of IRDC 18223-3. The inte- 
gration is from 30-60 km s which includes the entire veloc- 
ity range of the line. Dotted contours of the 1.2 mm single dish 
continuum map (Beuthe r et al. 2002at are overlaid and shaded in 
grayscale. The star marks the position of the 1.3 millimeter con- 
tinuum peak which is approximately 3.5' south of IRAS 18223- 
1243. We attribute the difference in peak position between the 
1.2 mm single dish and 1.3 mm interferometer data to pointing 
error in the single dish data. The size of the SMA primary beam 
(54") is marked by the dashed circle. The beam size is shown in 
the bottom right corner. 



3 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain) 



vations (Beut her et al. 20 02a) highlights the filamentary struc- 
ture, and the 12 CO(2-l) mostly traces emission associated with 
the outflow. We also detect a weaker peak approximately 40" 
east and slightly north of the primary peak. However, as it is out- 
side the primary beam of the interferometer observations, we do 
not detect this peak in subsequent observations discussed here. 
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Fig. 2. Integrated flux of IRDC 18223-3 observed with PdBI of the 3mm dust continuum, SiO, HN 13 C, and H 13 CO + . The integrations 
include the entire velocity range of the lines, namely 35-54 km s~' for SiO, and 43-48 km s _1 for both HN 13 C and H 13 CO + . Contours 
always start at the 3cr level. For the continuum, SiO, HN 13 C, and H 13 CO + respectively, contours increase in steps of 3, 5, 3, and 
8 cr where cr is 0.15, 11, 16, and 11 mJy birr 1 . The velocity (1st) and line width (2nd) moment maps of HN 13 C and H 13 CO + are 
presented in the panels on the right. The redshifted CO emission contours are overlaid on the H 13 CO + 2nd moment map. The star 
indicates the position of the dust continuum peak. The size of the synthesized beam is shown by the ellipse at the bottom left corner 
of each panel. 



It is interesting to note that this CO(2-l) peak is not associated 
with any emission in the 1.2 mm dust continuum. 

Figure [2] presents the PdBI observations of several molec- 
ular species along with the dust continuum. We defer discus- 
sion of SiO, a shock tracer, to section 13.31 Outflow Properties, 
because separate blueshifted and redshifted integrated maps 
(Figure |5j indicate that it is associated with outflow activity. 
HN 13 C and H 13 CO + both trace the large filamentary structure 
associated with IRAS 18223-1243 seen previously by Beuther 
et al. (2005c). Even the 3 mm dust continuum exhibits elonga- 
tion along the filamentary structure. The slight protrusion west- 
ward of the central 3 mm continuum peak may come about as 
a contribution from the outflow as there is also a peak in SiO 
at approximately the same spatial position. As we do not see 
any indication of large scale velocity gradients in the first (ve- 
locity) moment maps of HN 13 C and H 13 CO + , there likely is no 
significant movement happening along the filament although we 
cannot exclude the expansion of gas in the plane of the sky. In 
contrast to previous N 2 H + (l-0) data (Beuther et al. 2005b), we 
do not see an increase in line width directly centered on IRDC 
18223-3 in the second (line width) moment maps. Rather, we see 
increased line width to the northeast of the central peak. The po- 



sition of this peak in the second moment maps does not appear 
to be related to the outflow. 

3.2. Millimeter Continuum Emission 

The interferometer continuum data (Figure[3]l resolve the central 
structure into several peaks at a spatial resolution of ~4000 AU. 
The 1 . 1 mm dust continuum traces one main peak with 65 mJy 
birr 1 and three secondary peaks. A secondary peak of 17 mJy 
bnT 1 occurs 4.0" to the northeast (PA. =26°) and another peak 
of 22 mJy birr 1 is situated 3.2" to the southwest (PA. =261°) 
These secondary peaks as well as a 27 mJy birr 1 peak that is 
more removed from the main peak (7.9" with PA. =209°) are 
also detected at 1.3 mm, but an extension to the northwest in the 
1.1 mm continuum is not seen at 1.3 mm. While the secondary 
peaks may be separate sub-sources, the line data presented below 
indicate that the two near-by peaks to the northeast and south- 
west are likely part of a larger scale infalling and rotating struc- 
ture. The spatial resolution of the 3 mm data is worse, but even 
in the zoomed in continuum image of Figure [3] the extension 
tracing the filament toward the south is already evident. 

With the flux measured within the outer contour of the 1 .3 
mm dust continuum, we estimate the gas mass M and beam av- 
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eraged column density Nh-, of the region following the methods 
of Hildebrand (1983) and Beuther et al. (2002a, 2005b) . These 
estimates are under the assumptions that the dust continuum is 
optically thin at this wavelength, and that the gas to dust mass ra- 
tio is 100. For these calculations we assume a temperature of 20 
K and a dust opacity index /3 of 2 which corresponds to an opac- 
ity of 0.35 cm 2 g _1 at 1.3 mm. Although Beuther & Steinacker 
(2007) determine components at two different temperatures, the 
20 K we assume here is closer to the temperature of the compo- 
nent with the majority of the mass in IRDC 18223-3. We mea- 
sure a total integrated flux of 130 mJy and a peak flux of 41 mJy 
bm -1 . From these values we derive a total mass of 120 M and 
a beam averaged column density of 3.8 x 10 24 crrT 2 correspond- 
ing to a peak visual extinction, A v (mag), of ~4000, where A v — 
Ntf 2 /0.94 x 10 21 cnr 2 mag -1 ( [Frerking et al. 1982| l. Taking the 
uncertainties of the temperature and dust opacity into account, 
we estimate that the calculated mass and beam averaged column 
density for the compact structures detected here are accurate to 
within a factor of approximately 5. 

It must also be noted that the mass estimate is heavily af- 
fected by missing flux in the interferometer data. Comparing the 
flux from the 250 GHz single dish bolometer map (dotted con- 
tours in Figure [1} with the 230 GHz SMA continuum flux, we 
estimate that without short spacing information, the interferom- 
eter filters out about 80% of the flux. Hence we recover only 
about 20% of the mass in our estimate. As a detection limit, we 
calculate a mass of 5.2 M Q and a column density of 5.3 x 10 23 
cm -2 at the 3<x (5.7 mJy brrT 1 ) flux level. The interferometer 
data are not sensitive to structures less massive than this. 



3.3. Outflow Properties 

Figure [4] demonstrates that there is a large scale outflow in 
the northwest-southeast orientation. The morphology of the 
blueshifted CO component to the southeast is cone-shaped and 
likely is tracing the outer walls of the outflow cavity. There 
is also large, not very collimated blueshifted emission to the 
northwest of the source seen predominantly in the 30 m data. 
Although much less prevalent, we detect some redshifted emis- 
sion on both sides of the primary 1.3 mm continuum peak. The 
fact that we detect red- and blueshifted emission on both sides 
of the driving source indicates that the outflow is roughly in the 
plane of the sky because as the outflow cavity widens, the near 
edge would appear blueshifted while the far edge would appear 
redshifted. 

We detect SiO toward IRDC 18223-3 with the PdBI (See 
Figures [2] and |5). This tracer of shocked material, especially in 
molecular outflows (ISchilke et al. 199Tb , is consistent with the 
outflow orientation we determine from I2 CO(2-l) observations, 
and the southeastern SiO component coincides quite well with 
the southern edge of the blueshifted 12 CO cone-shaped emission. 
The slight difference in alignment between blue- and redshifted 
SiO emission shown in Figure [5] and the CO outflow orientation 
indicated by the arrows overlaid in this figure could potentially 
arise from the presence of multiple outflows. Such a scenario 
would indicate the presence of an undetected secondary source. 
We discuss the possibility of multiplicity further in section 13.61 
Alternatively, the difference in morphology between SiO and CO 
may be explained by precession of the jet/outflow in which the 
SiO emission traces the current jet component's location. The 
fact that SiO tends to trace the more collimated jet component 
whereas CO traces the less collimated outflow cavity provides 
support that this explanation is plausible. The CO data indicate 
that the outflow is in the plane of the sky, but a slight inclination 



is necessary to explain the fact that we detect blueshifted SiO 
emission predominantly to the northwest and redshifted emis- 
sion predominantly to the southeast indicating that the jet com- 
ponent is not exactly in the plane of the sky. In either case, the 
SiO data are consistent with the CO outflow to the extent that any 
disk component, if present, should have a position angle roughly 
45° east of north. 

Applying the methods of Cabrit & Bertout (1990) we de- 
rive properties of the system such as outflow mass, dynami- 
cal age, and outflow energy. These calculations assume that the 
I3 CO (2-l)/ 12 CO (2-l)line wing ratio remains constant over the 
entire outflow. Based on our single dish data which is consistent 
with the data of Choi et al. (1993), we adopt a value for this ratio 
of 0.1. Based on measuring a CO emission spatial extent of 25" 
from the central peak and adopting maximum velocities of 10.5 
km s" 1 and 13 km s _1 for the red- and blueshifted lobes respec- 
tively, the derived characteristics of the outflow associated with 
IRDC 18223-3 are listed in Tableg] 

The true velocity and spatial extent of an outflow must in- 
clude a correction for the inclination of the system from the 
plane of the sky. Thus any outflow parameters derived from the 
measured velocity and outflow extent also depend on this cor- 
rection. Based on the observation of blueshifted 12 CO emission 
on both sides of the central source, we interpret this to mean that 
the outflow is roughly in the plane of the sky. In this interpreta- 
tion, our measurement of the CO spatial extent is a reasonable 
approximation of the true outflow extent and the effect of incli- 
nation on the calculated outflow properties should be small. The 
velocity correction, on the other hand, may be significant for in- 
clinations close to 0°, leading to a much shorter timescale than 
that listed in Table [3] In turn, the outflow rate, mechanical force 
and mechanical luminosity may be significantly larger than the 
values given in Table[3]since these values depend on the outflow 
timescale. 

We compare the quantities presented in Table [3] with the sin- 
gle dish study of massive outflows in Beuther et al. (2002b). Our 
calculated mechanical force (F m ), outflow mass (M t ), outflow 
rate (M out ) and luminosity (L m ) place IRDC 18223-3 well within 
the massive regime of Figures 4a (F m vs. M C0le ) and 7 (M t vs. 
M core ) and at the lower end of the high-mass regime in Figures 4b 
(F m vs. L m ) and 5 (M out vs. L m ) of Beuther et al. (2002b). In all of 
these figures, IRDC 18223-3 is well above the low mass sources 
from the studies by Cabrit & Bertout (1992) and Bontemps et 
al. (1996). Both of the figures in which IRDC 18223-3 is situ- 
ated near the lower end of the high-mass regime are luminosity 
dependent. The position of IRDC 18223-3 in these plots is a re- 
flection of the source's youth and we expect that the luminosity 
will increase with time eventually shifting the source well into 
the regions occupied by the high-mass outflows in Figures 4b 
and 5. This general agreement with the properties of the high- 
mass sources further supports the interpretation that this out- 
flow is associated with a source in the massive regime. The 
outflow parameters derived here confirm previous suggestions 
(IBeuther et al. 2005bllBeuther & Steinacker 20071 1 that this truly 
is a massive star in the making. 

It is also interesting to note that in the SiO integrated map 
in Figure [2] there is an emission peak approximately 7" west 
of the main peak and Figure [5] demonstrates that this includes 
both red- and blueshifted velocities. We also detect a collimated 
redshifted CO feature situated westward of the primary dust 
continuum peak. Puzzlingly, we do not detect a corresponding 
blueshifted counterpart, nor a millimeter source associated with 
this CO feature, but it may have a connection to the SiO detec- 
tion as it appears to originate close to where SiO has a secondary 
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Fig. 3. The 1.1 mm (SMA), 1.3 mm (SMA), and 3.4 mm (PdBI) dust continuum. The area displayed for all three panels is identical. 
Contours start at 3<x and increase in steps of 3<x where cr is 2.2, 1.6, and 0.15 mJy birr 1 for the 1.1, 1.3, and 3.4 mm continuum 
images respectively. 



Table 3. Outflow parameters without inclination angle correction and corrected for an inclination of 20° with respect to the plane of 
the sky. 



if] 


M, [M e ] 


p [M km/s] 


E [erg] 


size [pc] 


t[yr] 


M out [Mo/yr] 


F m [M /km/s/yr] 


Lm [L©] 




13 


150 


1.7 x 10 4t> 


0.45 


37000 


3.5 x 10 4 


4.0 x 10~ J 


3.8 


20 


13 


440 


1.5 x 10 47 


0.48 


14000 


9.6 x 10~ 4 


3.2 x 10~ 2 


89 



Entries include inclination i, total outflow mass M,, momentum p, energy E, size, outflow dynamical age t, outflow rate M out , mechanical force F, 
and mechanical luminosity L m . Inputs are discussed in the text. 



peak. However, this is quite speculative, and we refrain from fur- 
ther comment and interpretation of this collimated redshifted CO 
emission. 

3.4. Age of the System 

Compared to other massive star formation regions, IRDC 18223- 
3 is extremely young as evidenced by the molecular proper- 
ties described in Beuther et al. (2005c). The C I8 line width 
measured at the main continuum peak of the SMA data is 1 .7 
km s~' whereas it is typically a factor of two to five times 
larger for molecules that trace similar densities in more evolved 
massive star formation regions (see e.g. [Bronfman et al. 19961 
IHatchell et al. 19981 IHofner et al. 20001 IBeuther et al. 2002al l. 
In the two secondary continuum peaks, the C ls O line width is 
smaller still. The C ls O line width at the peak to the northeast is 
0.94 km s _1 and in the southwest peak it is 0.60 km s -1 . While 
the line width at the center is less than 2 km s~', it is still larger 
than the line widths in the secondary peaks confirming that the 
outflow seen in 12 CO(2-l ) originates from the primary peak. The 
quite narrow central line width that we measure, along with the 
paucity of chemical species present toward the center, support 
the interpretation that IRDC 18223-3 is at a very young evolu- 
tionary stage. 

When the luminosity published by Beuther & Steinacker 
(2007) is considered in the framework of the models by 
Krumholz et al. (2007), estimates for the age of the system are 
less than ~10,000 years. Here we estimate the outflow dynami- 
cal age of IRDC 18223-3 to be approximately 37,000 years. 

This timescale calculation depends on the tangent of the out- 
flow's inclination angle (with 0° being in the plane of the sky). 




10 -10 -20 

AR.A. ["] 



Fig. 5. The SiO(2-l) integrated map. Again the triangle sym- 
bol marks the position of the dust continuum peak. The arrows 
indicate the orientation of the 12 CO(2-l) outflow as shown in 
Figure |U The PdBI beam is included in the lower right corner. 
Blueshifted emission (solid contours) integrated over velocities 
of 35-39 km s in steps of 4cr starting at 3cr where cr=0.011 
Jy birr 1 . Redshifted emission (dashed contours) integrated over 
velocities of 50-54 km s _1 in steps of 4cr starting at 3cr where 
<x=0.011 Jybm 1 . 



Applying a correction for inclination could potentially result in 
a much younger age than the 37,000 years which we calculate 



C. Fallscheer et al.: Rotational Structure and Outflow in IRDC 18223-3 7 




20 10 -10 -20 20 10 -10 -20 20 10 -10 -20 
AR.A. ["] AR.A. ["] AR.A. ["] 

Fig. 4. 12 CO(2-l) integrated maps of IRDC 18223-3 from SMA, IRAM 30 m, and SMA + IRAM 30 m combined in the left, 
center, and right panels respectively. Panels that include SMA data contain both extended and compact configuration observations. 
We define the outflow orientation (indicated by the arrows) by the blueshifted interferometer component. In each panel, the 1.3 
mm continuum peak is indicated by the triangle, and the beam is included in the lower right corner. Left: Blueshifted emission 
(solid contours) integrated over velocities of 38.5-43.5 km s in steps of 3cr starting at 3<x where cr=0.12 Jy birr 1 . Redshifted 
emission (dashed contours) integrated over velocities of 50-54 km s _I in steps of 2cr starting at 3cr where cr=0. 17 Jy brcT 1 . Center: 
Blueshifted emission (solid contours) integrated over velocities of 36.5-43.5 km s~' in steps of lcr starting at lcr where <x=6.4 Jy 
bm . Redshifted emission (dashed contours) integrated over velocities of 46.5-53.5 km s in steps of lcr starting at lcr where 
<x=6.9 Jy birr 1 . Right: Blueshifted emission (solid contours) integrated over velocities of 36-44 km s~' in steps of 6cr starting at 
3cr where cr-Q. 14 Jy birr 1 . Redshifted emission (dashed contours) integrated over velocities of 50-56 km s _I in steps of 6cr starting 
at 3cr where cr=0.14 Jy birr 1 . 



here. As seen in Table [3] for example, with an inclination angle 
of 20° with respect to the plane of the sky, the age estimate de- 
creases by a factor of nearly 3. 

Both the age estimate based on luminosity and the outflow 
dynamical age are subject to large uncertainties making it diffi- 
cult to pinpoint the true age. We are rather in support of interpret- 
ing the source simply as very young based on the non-detection 
at infrared wavelengths as well as its not very evolved chemical 
state. 

3.5. Rotational Structure 

In this section we present maps of the velocity distribution 
weighted by intensity of the molecules in our SMA data set that 
are associated with tracing the physical conditions in dense gas. 

The velocity moment map of C I8 (Figure[6]) exhibits a ve- 
locity gradient of approximately 2 km s _1 from south to north. 
The observed gradient is neither aligned with nor orthogonal to 
the outflow orientation, presumably because this tracer of denser 
regions is being strongly affected by infall and the outflow itself. 
The position angle of the observed gradient differs by about 30° 
with respect to the 12 CO outflow orientation. As C 18 is an iso- 
topologue of 12 CO, it may indeed be tracing a disk-like structure 
although the observed gradient may be influenced by infall and 
the outflow. This is not an entirely unexpected result despite the 
fact that C ls O is optically thin and not a good outflow tracer. 

Methanol is well documented as a molecule that 
traces cores, shocks, and masers in star formation re- 
gions (e.g., |j0rgensen et al. 2004{ IBeuther et al. 2 005a; 



Sobolev et al. 2007]). It has also been reported as a low 
mass disk tracer (e.g. [Goldsmith et al. 19991) . yet it has not 



previously been associated with tracing disks in massive star 
formation. In the velocity moment map of the CH30H(6o-5o) 
line at 290.110 GHz (middle panel of Figure |6]), we see a 
flattened, elongated rotating structure perpendicular to the 
outflow axis. The structure evident here is consistent with the 
velocity moment maps of the next two strongest methanol lines, 
229.759 GHz and 290.070 GHz. The remaining methanol lines 
in our spectral setup (listed in Table [TJ are too weak to be able 
to justify similar statements. The approximately 7.5" extent 
corresponds to a major axis extent of 28,000 AU at the given 
distance of 3.7 kpc. Although this is rather large to be a disk, the 
narrow velocity range of 3 km s _1 coupled with the large spatial 
extent perpendicular to the outflow support an interpretation that 
this is a large rotating and infalling core similar to the toroids 
described in the recent review by Cesaroni et al. (2007). There 
is also a high velocity component in the northwest corner of the 
CH3OH velocity moment map which coincides with the outflow 
and is likely associated with the extension seen in the 1 . 1 mm 
dust continuum overlay of Figure [6] 

In Figure|7] we present the position-velocity diagrams of the 
CH3OH and N2H + emission. The cuts go through the peak of the 
dust continuum and have position angles along the velocity gra- 
dients seen in the moment maps. In the N2H + position-velocity 
diagram, it is apparent that we suffer from missing flux espe- 
cially around the velocity of rest and at more redshifted veloci- 
ties, but we still see an indication for the same trends as in the 
CH3OH position-velocity diagram. Because of the larger effect 
of missing flux in N2H + compared to CH3OH (see Figure |7), 
we do not further analyze this molecule suited for probing early 
stages of star formation, but rather concentrate on CH3OH. 



8 



C. Fallscheer et al.: Rotational Structure and Outflow in IRDC 18223-3 



Although the position-velocity diagrams do not exhibit 
Keplerian velocity signatures, they are consistent with the veloc- 
ity moment maps in the sense that we see the largest velocities 
furthest from the center. In Keplerian motion, the equilibrium 
condition between gravitational and centrifugal forces means 
that the velocities at large distances from the source should be 
lower than those closer to the center. One explanation for the de- 
viation from Keplerian motion is that we do not see all the way 
to the center due to optically thick emission, and hence see only 
rotation of the outer envelope. However, without optical depth 
information, we cannot assume this to be the case, so we in- 
terpret the observed velocity gradient simply as non-Keplerian 
rotation of a large structure perpendicular to the outflow orienta- 
tion. This should not come as a surprise as the mass of the central 
source is likely much less than the mass of the rotating structure. 

3.6. Theoretical Modeling 

Assuming optically thin line emission, we modeled the velocity 
gradient seen in the CH3OH velocity moment map. The model- 
ing was done using the UMch Model dUlrich 19761 IKeto 2007t 
which is an angular momentum conserving mass inflow model 
that includes gravity and an isothermal rotating envelope. We 
calculated the velocity and density along the line of sight and 
then multiplied these together to produce the velocity distribu- 
tion. The model does not include a disk because this adds an- 
other parameter which we cannot distinguish between due to 
observational limitations. The input parameters we used were 
determined from the CH3OH moment map. We measure an ob- 
served rotation speed of 1.5 km s _1 at 14,000 AU (3.8" at the 
given distance of 3.7kpc). We chose a mass infall rate of 10~ 3 M 
yr _1 . The system in our current studies does not fit well into the 
picture of Keplerian rotation that describes low-mass accretion 
disks (Dutrey et al. 2007 ), but in order to get a rough estimate of 
the mass in the rotating structure we make the assumption that 
rotational and gravitational forces are balanced at the outer edge 
of the disk. We then calculate the mass that is in gravitationally 
bound rotational motion within a given radius, r, using the fol- 
lowing equation: 



M rof [Mo] = 1.13 ■ 10~ 3 x <5v 2 [km/s] x r[AU] (2) 

We take r to be 14,000 AU, and 6v to be 1.5 km s _1 , half the ex- 
tent of emission and half the velocity range seen in the CH3OH 
velocity moment map. From these values, we estimate M ro , to 
be 36 M Q . The rather low luminosity of this source leads to the 
interpretation that this mass is strongly dominated by a larger ro- 
tational structure opposed to a central object. In addition to the 
large uncertainties (typically within a factor of 5) associated with 
mass estimates in general, the discrepancy between the mass es- 
timate of 120 M from the continuum emission (Section |3.2t 
and the mass calculated here can be explained by the fact that 
the system is not rotationally supported but rather infalling, as 
we assume for the modeling. 

The resulting density distribution and velocity field along 
with the combination of these into the velocity distribution are 
shown in Figure [8] In these plots the disk formation radius, Rq, 
is the radius in the midplane that has the highest density and cor- 
responds to the point at which the rotational force balances the 
gravitational force. At the disk formation radius, the Keplerian 
velocity is 2.1 km s _1 . We do not observe Keplerian velocity 
signatures in the position-velocity diagrams (Figure[7]l, so we do 



not expect to measure the Keplerian velocity predicted by the 
model. Although we do not exactly fit our quantitative param- 
eters, this simple model of rotationally supported infall repro- 
duces our observations well. 

In order to make a more direct comparison between the the- 
oretical model and our observations, we have convolved the 
model to the 1.7" angular resolution of the data and rotated it 
counterclockwise by 45° such that it has the same orientation as 
the observed velocity gradient. Figure [9] presents the smoothed 
and rotated model and the CH3OH velocity moment map side by 
side. We emphasize that this model is not meant to be a best fit 
of the data, but rather to show that our interpretation presented 
here is physically plausible. 

Keplerian rotation signatures are not detected in 
this and other regions of massive star formation (see 
ICeseroni et al. 2007j l. One scenario that could explain this 
is that as the core collapses, it fragments into smaller con- 
densations and forms several stars, each with accretion disks 
not necessarily aligned with the parental core's rotational 
axis. If this were the case, we would expect to see a more 
chaotic velocity structure as well as more outflow signatures. 
Additionally, this may not be physically reasonable in light of 
the modeling by Krumholz et al. (2007, 2009). In these models, 
multiple stars form within the same accretion disk meaning that 
the rotational axis should be similar for all sub-sources. Since a 
driving source is required to power the outflow, we suspect that 
a true Keplerian accretion disk may have formed at the center 
on scales we cannot observationally resolve. 

One interpretation of the observed lobes to the northeast and 
southwest in the millimeter continuum is that there are several 
unresolved sources in the vicinity. If these objects are all part of 
an inward spiraling structure, this might be mistaken for a veloc- 
ity gradient across them. However, our ability to reproduce the 
gradient in the CH3OH velocity moment map through successful 
modeling provides evidence against this interpretation and rather 
in favor of a scenario in which a single component is surrounded 
by a large rotating object. In order to discriminate between these 
scenarios mentioned here, even higher spatial resolution images 
of the system would be necessary. Further studies of spectral 
lines probing deeper into the core would also likely prove to be 
helpful. 



4. Conclusions 

The Infrared Dark Cloud 18223-3, a source at the onset of mas- 
sive star formation, has been observed with multiple submillime- 
ter observatories. Despite being at a very early stage in the evo- 
lutionary process, we detect a molecular outflow and evidence 
for a large rotating object perpendicular to the outflow. 

This region likely represents one of the earliest stages of disk 
formation. IRDC 18223-3 exhibits a large rotating and flattened 
envelope structure which we suspect is one of the earliest de- 
tectable stages in the evolutionary sequence of accretion disks. In 
this framework, the structure we see and have described here will 
continue contracting into a true accretion disk as time passes. 

The combined SMA and 30 m data reveal a well defined 
northwest-southeast outflow orientation. To the southeast of the 
1.3 mm continuum peak, we see a cone-like blueshifted com- 
ponent and the northwest is dominated by a broader blueshifted 
component. While it is less clear what is happening in the red- 
shifted regime, the blueshifted components indicate that we are 
likely looking at a system that is roughly in the plane of the sky. 
From this we deduce that the likely disk orientation would be 
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Fig. 6. Upper Row: C ls O(2-l), CH 3 OH(6 -5 ) A, and N 2 H + (3-2) total intensity (0th) moment maps. Lower row:C ls O(2-l), 
CH30H(6o-5o) A, and N2H + (3-2) velocity (1st) moment maps. The high velocity component toward the northwest corner of the 
CH3OH velocity moment map is also seen as an extension in the 1 . 1 mm dust continuum overlaid on the middle panel. The contours 
start at 5<x and increase in steps of 5cr where cr is 2.25 mJy birr 1 . All moment maps were clipped at the six sigma level of the 
respective line's intensity map. The star symbol indicates the position of the 1.3 mm dust continuum peak. The beam size is shown 
in the lower left corner of each plot. 
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Fig. 7. Position-velocity diagrams derived for cuts along the observed velocity gradient. The offset refers to the distance along the 
cut from the dust continuum peak. The rest velocity of 45.3 km s _1 is marked in both plots. Contours for both diagrams start at 0. 16 
Jy and increase in 0. 16 Jy steps. 



10 C. Fallscheer et al.: Rotational Structure and Outflow in IRDC 18223-3 



kms 

44 45 46 47 44 45 46 47 




4 2 -2 -4 4 2 -2 -4 

AR.A. ["] AR.A. ["] 



Fig. 9. Left: The model shown in the lower panel of Figure [8] after convolving to the 1.7" resolution of the data and rotating 
counterclockwise by 45° into alignment with the observed velocity gradient. Right: The CH3OH velocity moment map as in Figure 
[6] The velocity range has been modified to exclude the high velocity components seen in the northwest of Figure [6] so that a more 
direct comparison to the model can be made. 



approximately edge on and that the associated velocity gradient 
would be in the northeast-southwest direction. 

Based on the observed outflow properties, we find that cal- 
culated parameters are consistent with a massive driving source. 
Although the central source is not very massive at this point in 
time, the outflow characteristics provide a strong indication that 
this star will evolve into the massive regime. 

Indeed we do see velocity gradients in CH3OH and N2H + , 
although they are observed over a very large spatial distance 
and a relatively small velocity range. On the order of 20,000 
AU in size, the large rotating core we currently detect is much 
larger than other disk or disk-like structures around similar low 
luminosity intermediate to high mass star forming regions pub- 
lished in the literature. For instance, the disk associated with 
IRAS 20126+4104 is less than a t hird the size of the struc- 
ture we observe in IRDC 18223-3 (Cesaroni et al. 2005 ). The 
reported disks aroun d G192.16-3.82 (|Shepherd et al. 2001| l and 
IRAS 18089-1732 dBeuther et al. 2004t are smaller still We 
thus suspect that CH3OH and N2H + are actually tracing the outer 
edges of an infalling toroid. This rotating toroid likely plays an 
important role in feeding an accretion disk within the unresolved 
central region and may decrease in size over time as the outer 
edges dissipate or contract leading to growth of the disk. Hence 
we see here potentially the earliest stages of the disk formation 
process. As there is not much known about CH3OH as a tracer of 
massive disks, one possibility is that it may be a good tracer for 
disk kinematics only at very early evolutionary stages, but that 
other hot core molecules such as CH3CN, for example, may be 
better for more evolved massive star formation regions. 

Using the Ulrich Model for mass infall we have successfully 
reproduced the velocity gradient observed in methanol. This 
agreement between model and observations further enhances the 
argument that the rotation we see results from a single flattened 



rotating entity rather than multiple sources or other complicated 
scenarios. 

Similarly, we detect a velocity gradient in C I8 that is not 
aligned with the 12 CO outflow axis. However, the spatial scale is 
large, and the velocity range is even smaller than in CH3OH and 
N2H + . In this case, since C ls O is an isotopologue of the standard 
outflow tracer, I2 CO, it is likely that the velocity gradient we see 
is being influenced by infall motion or the outflow itself which 
would tend to shift the observed gradient toward the direction of 
the outflow orientation. 

The results we have presented further support the idea that 
IRDCs are sources undergoing the early stages of high-mass star 
formation. Thus IRDCs provide ideal sites to study the birth of 
massive stars. 

Acknowledgements. C.F. and H.B. acknowledge support by the Deutsche 
Forschungsgemeinschaft, DFG project number BE 2578. C.F. also acknowledges 
support from the International Max-Planck Research School for Astronomy &■ 
Cosmic Physics at the University of Heidelberg and thanks Jan Martin Winters 
for his guidance in reducing the PdBI data sets. We express our gratitude to the 
anonymous referee for his/her constructive comments which improved the paper. 

References 

Beuther, H., Schilke, P., Menten, K. M., et al. 2002a, ApJ, 566, 945 

Beuther, H., Schilke, P., Sridharan, T.K., et al. 2002b, A&A, 383, 892 

Beuther, H., Hunter, T. R., Zhang, Q., et al. 2004, ApJ, 616, L39 

Beuther, H., Zhang, Q., Greenhill, L. J., et al. 2005a, ApJ, 632, 355 

Beuther, H., Schilke, P., Menten, K. M., et al. 2005b, ApJ, 633, 535 

Beuther, H., Sridharan, T.K., & Saito, M. 2005c, ApJ, 634, L185 

Beuther, H., Churchwell, E. B., McKee, C. F, & Tan, J. C. 2007, in Protostars 

and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. Arizona 

Press), 165 

Beuther, H., & Steinacker, J. 2007, ApJ, 656L, 85 
Beuther, H., & Walsh, A. J. 2008, ApJ, 673L, 55 

Bontemps, S., Andre, P., Terebey, S., & Cabrit, S. 1996, A&A, 311, 858 



C. Fallscheer et al.: Rotational Structure and Outflow in IRDC 18223-3 



11 




Hofner, P., Wyrowski, R, Walmsley, C. M., & Churchwell, E. 2000, ApJ, 536, 
393 

J0rgensen, J. K., Hogerheijde, M. R., van Dishoeck, E. E, et al. 2004, A&A, 415, 
1021 

Keto, E. 2007, ApJ, 666, 976 

Krumholz, M., Klein, R., & McKee, C. 2007, ApJ, 656, 959 
Krumholz, M., Klein, R., McKee, C, Offner, S., & Cunningham, A. 2009, 
Science, 323, 754 

Schilke, R, Walmsley, C. M., Pineau des Forets, G., & Flower, D. R. 1997, A&A, 
321,293 

Shepherd, D. S., & Churchwell, E. 1996, ApJ, 457, 267 
Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2001, Sci, 292, 1513 
Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2002, ASPC, 267, 415 
Sobolev, A. M., Cragg, D. M., Ellingsen, S. P., et al. 2007, IAUS, 242, 81 
Sridharan, T.K., Beuther, H., Saito, M., Wyrowski, F., & Schilke, P. 2005, ApJ, 

634, L57 
Ulrich, R. 1976, ApJ, 210, 377 
Yorke, H., & Sonnhalter, C. 2002, ApJ, 569, 846 
Zhang, Q., Hunter, T. R., Brand, J., et al. 2001, ApJ, 552, L167 
Zhang, Q., Hunter, T. R., Brand, J., et al. 2005, ApJ, 625, 864 
Zinnecker, H, & Yorke, H. 2007, ARA&A, 45, 481 



Fig. 8. Results of modeling the velocity gradient seen in CH3OH 
with the Ulrich infall model. Rj is the radius, and 1 Rj corre- 
sponds to 10,000 AU. Upper: The density distribution overlaid 
with the velocity field vectors. The velocity vectors vary between 
0.5 and 1.0 km s _1 in length, and the densities vary between 
9.4xl0 5 and 3.9xl0 7 cm 3 . The logarithmic plotting of the den- 
sity resembles a disk, but we have not included one in this mod- 
eling. Lower: The velocity and density distributions above mul- 
tiplied together then averaged along the line of sight to reproduce 
the velocity gradient seen in the CH3OH velocity moment map. 



Bronfman, L., Nyman, L.-A, & May, J. 1996, A&AS, 115, 81 
Cabrit, S., & Bertout, C. 1990, ApJ, 348, 530 
Cabrit, S., & Bertout, C. 1992, A&A, 261, 276 

Cesaroni, R., Felli, M., Testi, L., Walmsley, C. M., & Olmi, L. 1997, A&A, 325, 
725 

Cesaroni, R., Neri, R., Olmi, L., et al. 2005, A&A, 434, 1039 

Cesaroni, R., Galli, D., Lodato, G., Walmsley, C. M., & Zhang, Q. 2007, in 

Protostars & Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. 

of Arizona Press), 197 
Choi, M., Evans, N. J., & Jaffe, D. T. 1993, ApJ, 417, 624 
Dutrey, A., Guilloteau, S., & Ho, P. 2007, in Protostars & Planets V, ed. B. 

Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. of Arizona Press), 495 
Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590 
Goldsmith, P., Langer, W., & Velusamy, T. 1999, ApJ, 519, L173 
Hatchell, J., Thompson, M. A., Millar, T. J., & MacDonald, G. H. 1998, A&AS, 

133, 29 

Hildebrand, R. H. 1983, QJRAS, 24, 267 



